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Abstract 

We consider the effect of the internal nuclear quadrupole interaction on quantum tunneling in 
complex multi-atomic two-level systems. Two distinct regimes of strong and weak interactions 
are found. The regimes depend on the relationship between a characteristic energy of the nuclear 
quadrupole interaction A* and a bare tunneling coupling strength Ao- When Ao > A*, the internal 
interaction is negligible and tunneling remains coherent determined by Ao- When Ao < A*, 
coherent tunneling breaks down and an effective tunneling amplitude decreases by an exponentially 
small overlap factor rf <C 1 between internal ground states of left and right wells of a tunneling 
system. This affects thermal and kinetic properties of tunneling systems at low temperatures 
T < A*. The theory is applied for interpreting the anomalous behavior of the resonant dielectric 
susceptibility in amorphous solids at low temperatures T < 5mK where the nuclear quadrupole 
interaction breaks down coherent tunneling. We suggest the experiments with external magnetic 
fields to test our predictions and to clarify the internal structure of tunneling systems in amorphous 
solids. 

PACS numbers: 6143.Fs, 77.22. Ch,75.50Lk 
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I. INTRODUCTION 



A transition between two energy minima separated by a potential barrier U occurs in 
different ways at high and low temperatures. At high temperatures the motion between the 
energy minima is classical, i.e., thermally activated, suggesting that the particle acquires 
the energy U from the environment to overcome the barrier. This results in an exponential 
Arrhenius factor exp(— U/T) for the transition rate. The classical transition rate decreases 
rapidly with decreasing the temperature. At low temperature the above-described transport 
crosses over into the weakly temperature dependent quantum tunneling. A tunneling tran- 
sition rate is essentially governed by the exponentially small and temperature-independent 
factor exp(— 2S(U)/h), where S(U) describes the classical action for the motion through the 
inverted potential barrier. The strong sensitivity of the tunneling exponent against particle 
mass makes quantum tunneling more favorable for light particles as electrons, while the 
motion of heavy nuclei is classical down to temperatures of order of Kelvin. However, for 
T < IK quantum tunneling displays itself in the thermodynamic and kinetic properties of 
atomic systems. One impressive example is amorphous solids, in which the low temperature 
properties are governed by the two-level systems (TLS's). These TLS's are made of atoms 
or groups of atoms, experiencing a tunneling motion between pairs of energy minima sepa- 
rated by potential barriers [1]. They contribute to the universal thermodynamic and kinetic 
properties in all known glasses as well as some other disordered materials for T < IK (see 
Refs. [2-4]). 

A TLS is described by the standard pseudospin 1/2 Hamiltonian 

h TLS = -A • s x - A • s z . (1) 

Here Ao is a tunneling amplitude coupling two energy minima and A is a level asymmetry. 
The quantity A is exponentially sensitive to external and internal parameters of the system. 
This results in a logarithmically uniform distribution of tunneling amplitudes for various 
TLS's 

P(A,A ) = -£-, P = const. (2) 

This distribution of TLS's leads to the universal temperature and time dependencies for 
various physical characteristics of amorphous solids, thus making "glassy" behavior very 
easily recognizable. For instance, this includes the logarithmic relaxation in time of the spe- 
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cific heat [1-4] and the non-equilibrium dielectric susceptibility [5]. In addition, dielectric 
and acoustic (a sound velocity) susceptibilities in glasses show a logarithmic temperature 
dependence also associated with the distribution Eq. (2). Note that the exponential sensi- 
tivity of a tunneling strength A to environmental interactions also broadens noticeably the 
distribution of tunneling amplitudes and relaxation rates in other systems, e.g., tunneling 
of a large electronic spin in magnetic molecule Mn^Ac. The latter system shows a broad 
spectrum of relaxation times due to interaction of electronic spins with nuclear spins [6, 7]. 
A broad spectrum of relaxation times also has been reported in new disordered magnetic 
alloys [8]. 

The nature of the tunneling systems in amorphous solids remains unclear in spite of 
the theoretical efforts attempting various models [9-13]. The main problem in the theory 
is a lack of experiments which are capable of proving the advantage of a specific model 
against the original phenomenological model [1] based on Eq. (2). The phenomenological 
approach Eq. (2) can be used to explain a variety of experimental data, while the interaction 
between TLS can successfully be treated as a weak correction [13, 14]. Even the internal 
TLS structure is unclear yet. For instance, nobody knows how many atoms do participate 
in a single tunneling event. We hope that understanding the nuclear quadrupole interaction 
effects in glasses, which are sensitive to the internal TLS structure, will help to resolve this 
question. 

Recent experimental investigations of amorphous solids at very low temperatures have 
revealed a number of qualitative deviations from the predictions of the standard tunneling 
model Eq. (2). In particular, it is demonstrated in several works [13, 15, 16] that, for 
T < 5mK, the expected logarithmic temperature dependence of the dielectric constant 
breaks down and the dielectric constant becomes approximately temperature- independent. 
This result conflicts with the logarithmically uniform distribution Eq. (2) of TLS's over 
their tunneling amplitudes. To resolve the problem, one can assume that the distribution 
of TLS's has a low-energy cutoff at A , min ~ 5mK. This assumption, however, contradicts 
the observation of very long relaxation times in all glasses. These times (a week or longer) 
require much smaller tunneling amplitudes [5] than 5mK (remember that the TLS relaxation 
time is inversely proportional to its squared tunneling amplitude). 

We suggest the explanation of this controversial fact by using the recent model of Wiirger, 
Fleischmann, and Enss [17] who proposed that the nuclear quadrupole interaction affects the 
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properties of TLS's at low temperatures by the mismatch of the nuclear quadrupole states 
in different potential wells (see Fig. 1). This is very similar to the electronic spin tunneling 
suppression by the nuclear spin interaction in magnetic molecules [7]. The significance of the 
nuclear quadrupole interaction has recently been proven experimentally in glycerol glass [18]. 
This interaction helps to understand the anomalous magnetic field dependence of dielectric 
properties in entirely non-magnetic dielectric glasses [16, 19, 20]. 

We show that an effective tunneling amplitude of TLS having an energy less than its 
nuclear quadrupole interaction is remarkably reduced due to the mismatch of TLS nuclear 
quadrupole states in its two energy minima, similarly to the polaron effect. Consequently 
the spectrum of TLS tunneling amplitudes A possesses a pseudogap below the nuclear 
quadrupole interaction energy A*. Dielectric and acoustic susceptibilities of glasses become 
temperature indepedent for temperatures belonging to this pseudogap because there is no 
TLS with A ~ T to contribute. Thus our model explains the experimental observations. In 
addition, we predict that the application of a strong external magnetic field will reduce the 
mismatch of different nuclear quadrupole states and thus restore the logarithmic temperature 
dependence of TLS susceptibilities. According to our theory dielectric glasses having no 
nuclear quadrupole interaction should obey the predictions of the tunneling model, which 
mostly agrees with the experiment. 

Since a low temperature dielectric constant and a speed of sound in glasses have a similar 
physical nature they should have a similar behavior at low temperatures. Therefore the 
saturation in the logarithmic temperature dependence of a sound velocity should be seen 
at T < lOmK in materials possessing non- vanishing nuclear quadrupole moment. The 
situation with sound velocity measurements is, however, more complicated because it is 
more difficult to perform the low-temperature measurements for a sound velocity then for 
a dielectric constant. There exists a few measurements that we discuss together with the 
measurements of the dielectric constant. The present theory is developed for the dielectric 
constant, although our conclusions can be extended to the velocity of sound without major 
changes. 

The paper is organized as follows. In Section II the model of two-level systems affected 
by the nuclear quadrupole interaction is introduced. Then, we discuss qualitatively the 
renormalization of tunneling amplitude by the nuclear quadrupole interaction depending on 
the relation between tinneling splitting A and a quadrupole interaction A* . In Section III 



4 



the expression for the resonant dielectric susceptibility is obtained and the influence on the 
tunneling amplitude renormalization on the resonant part of the TLS dielectric constant is 
described in the qualitative level. In Section IV we use a perturbation theory to characterize 
quantitatively the temperature dependence of the dielectric constant in the presence of 
nuclear quadrupole interactions in the high and low temperature limits. In Section V we 
consider a solvable model of TLS's coupled to harmonic oscillators, conveniently replacing 
nuclear spins. The solution of this problem permits us to simulate the influence of the 
nuclear quadrupole interaction on a TLS dielectric constant at all temperatures of interest. 
In Section VI the effect of the external magnetic field on the dielectric constant is considered 
within the simplified model of Sect. V. In Section VII the parameters of our model are 
compared with the experimental data. In final Section VIII the conclusions are formulated 
and the suggestions for an experimental verification of our theory are made. The short 
version of the manuscript appears in the Physical Review Letters [21]. 

II. MODEL 

A. Nuclear Quadrupole Interaction 

How does the nuclear quadrupole interaction affect tunneling? Consider a tunneling 
system formed by n atoms all possessing a nuclear spin I > 1 and consequently a nuclear 
electrical quadrupole moment. The total tunneling Hamiltonian h can be described by the 
standard TLS pseudospin Hamiltonian Iitls (1) an d the quadrupole interactions H r in the 



The local nuclear quadrupole Hamiltonians H r i can be expressed as a sum of interactions 
of all n nuclear spins Ii over all n atoms that simultaneously participate in its tunneling 
motion with the local electric field gradient tensors dFa' r ^ /dxb different in general for the 



right well (s z = 1/2) and Hi in the left well (s 2 



1/2) as follows 




(3) 
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right and left wells 

n 



i=l 

a,b=x,y,z 

Here if is a nuclear spin projection onto the a— axis and Q is the electrical quadrupole 
moment. 

In what follows, we consider a simplified model for the nuclear quadrupole interaction 
(4) possessing the axial symmetry (see Fig. 1) 



if,, r = f.((C) 



» /(/ + !)-, (5) 
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Here 

6wQ|9F/9a;|, (6) 

and ux, and define the directions of the electric field gradient in the left and right wells, 
respectively. 

z A u L r\ < z 




FIG. 1: The two level system having different nuclear quadrupole quantization axes ul and ur 
defined by the local electric field gradient in the left and right wells. 



There exists also the magnetic-dipolar interaction of nuclear spins. It is usually by several 
orders of the magnitude smaller than the quadrupole interaction and can be neglected. 

B. TLS ground state and its effective tunneling amplitude 



It will be shown below that the resonant dielectric susceptibility of a TLS ensemble is 
due to the ground state of tunneling systems affected by the external electrical field. Only 
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tunneling systems for which A ~ A can contribute to the dielectric constant. These two 
parameters must exceed the thermal energy, i.e., A > T (see [3, 4] and Sect. Ill for details). 
The structure of the ground state of the Hamiltonian (3) depends on the relation between 
typical tunneling amplitude A ~ T and characteristic value of the quadrupole interaction 
nb. 

If A > nb, the nuclear quadrupole interaction can be treated as a small perturbation 
and the TLS behavior obeys the standard tunneling model. 

In the opposite case A < nb, the tunneling term can be treated as a weak perturbation. 
In a zero-order perturbation theory the states in the left and right well can be considered 
separately. These states form a proper basis for eigenstates of a tunneling system. In partic- 
ular, the actual ground state is a linear superposition of the ground states in the isolated left 
and right wells. The energies corresponding to these ground states are E g i = —A/2 + E g ^, 
E gr = A/2+E g>r , where E g>r , E g j are the ground state energies of the quadrupole interaction 
Hamiltonians H r j (Eq.(3)) in the right or left wells, respectively. 

The effective tunneling amplitude A *, which couples any pair of levels in the left and 
right wells of TLS's with A < nb decreases due to a mismatch of nuclear quadrupole ground 
states in two potential wells. Obviously, when H r = Hi, the nuclear quadrupole interaction 
does not affect tunneling. This mismatch can be expressed in terms of a characteristic 
overlap integral between the nuclear spin ground states in the two wells. In the absence of 
tunneling one has 

v* = v n , v = (h I r g), (7) 

where the symbols \lg > (\rg >) stand for the ground-state wavefunction of the nuclear spin 
of a single tunneling atom in the left (right) well. This overlap integral enters directly into 
the effective tunneling amplitude. In fact, at very low temperatures and small tunneling 
amplitude A < nb one may treat the TLS in terms of the ground states in the right and 
left wells. The tunneling matrix element A * between these two states resulting from the 
perturbation — A s x is 

A * = n*A . (8) 

Thus, the overlap integral 77* determines the reduction of the tunneling amplitudes due to 
the nuclear quadrupole interaction. 

The specific value of the overlap integral depends on the absolute value of the nuclear 
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spin, nature of the nuclear quadrupole interaction, the number of tunneling atoms per one 
TLS, and the difference in the electric field gradients in the right and left wells. Currently, 
reliable information about all these parameters is not available. In principle the value of the 
overlap integral r) fluctuates from atom to atom. This will lead to the log-normal distribution 
of the overall overlap integral 77* in the large n limit. As we will see below (e. g. Eq. (56)) 
the behavior of the dielectric constant is mostly sensitive to In (77*), which possesses the 
gaussian distribution, so we can use its average value with the logarithmic accuracy of our 
consideration. 

To estimate typical value r/* let us turn to a simplified model for the nuclear quadrupole 
interaction Eqs. (4), (6), possessing the axial symmetry. Let be the angle between the 
directions of the electrical field gradient in the left and right wells. 

First, consider the case of an integer nuclear spin. If the interaction constant is positive, 
i.e., b > 0, the energy minimum is obtained when the projection = 0. In this case the 
ground state is non-degenerate. This simplifies our consideration and makes it possible to 
calculate the overlap integral of the ground state wavefunctions depending on the mismatch 
angle <fi (see Fig. 1). We confine ourselves with two simple cases of integer nuclear spins, 
namely I = 1,2. Then the overlap integral between the left- and right ground states can be 
expressed as 



We will pay most attention to these two non-degenerate cases because they are very conve- 
nient for the investigation of the effect of an external magnetic field on the overlap integral 
(see. Sect. VI.) 

The value of the rotational angle in glasses is unknown. We expect that it can be 
estimated by extended molecular dynamics simulations [22]. On the other hand, in the 
orientational glass (KBr)i_a; (KCN)^ (see Ref. [23]) the CN group rotates between different 
equilibrium positions by angle <fi = cos~ 1 (l/3) (see Ref. [3]). We assume that in glasses 
the rotational angle is the same and neglect its fluctuations due to a structural disorder. 
Then, for spin 1 = 1, one can estimate the overlap integral Eq. (9) as 



V = cos(0), 1=1; 
77 = |cos 2 (0) — sin 2 



2 (0)/2| ,1 = 2. 



(9) 



77 « cos(0) « 0.33. 



(10) 
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When the TLS contains n atoms tunneling simultaneously, the characteristic overlap inte- 
grals becomes 

r}* = (l\r) &r} n . (11) 

To understand and interpret experimental data, we assume that the total overlap integral 
is small 

?7* < 1. (12) 

This assumption is justified by the exponential decrease of the overlap integral with the 
number of atoms n participating in a single tunneling system (TLS). 

The case of a half-integer spin is more complicated because of the Kramers degeneracy. 
The quadrupole spin Hamiltonian of a tunneling atom i has two orthogonal ground states 
\il+ >, \il— > in the left well and two orthogonal ground states \ir+ >, \ir— > in the right 
well. As a result, the ground state of n tunneling atoms is 2 n -fold degenerate. In order to 
reduce the problem to a single pair of levels (one in the left well and the other in the right 
one), we divide the degenerate ground state into 2 n pairs of left \il > and right \ir > ground 
states coupled only with each other. In fact, tunneling amplitude between left and right 
states of a tunneling atom i is directly proportional to the overlap integral of the two states 
involved < il\ir >. Let us introduce two superposed states in the left well as follows 

| 11) = cos(a) | il+) + sin(a) | il—), 
| 12) = — sin(a) | il+) + cos(o;) | il—) 

The angle a is defined by the condition of orthogonality 

= (ir+ | 12) = (ir- \ 11), 
= —sm(a){ir+ \ il+) + cos(a)(ir+ \ il—), 
= cos(a) (ir— \ il+) + sin(a)(ir— \ il—). 

This equation can be solved for a if the equation determinant is equal to zero, i. e. 

| ir+)(il- | ir+) + \ ir-)(il- \ ir-) = 0. (15) 

This condition is satisfied when pairs of eigenstates \ir+ >, \ir— > and \il+ >, > are 
orthogonal to each other. One can verify the validity of Eq. (15) by projecting the left 
vectors \il+ >, \il— > onto the right vector subspace (\ir+ >, \ir— >). Choosing the basis 
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(13) 



(14) 



for the left well from the pairs of states Eq. (13), we can construct 2 n possible states of 
nuclear spins from the products of different left states (two choices for each atom). The basis 
states in the right well can be constructed similarly. Then each of the 2™ left basis states 
possesses a non-zero overlap integral with the only single single state from the 2 n right basis 
states. 

Thus, the degenerate states form groups of pairs of states connected by tunneling, while 
the states belonging to different pairs are uncoupled. Since all 2™ pairs possess identical 
properties with respect to tunneling, one can treat them as independent pairs of states. 
Then, our considerations become similar to those for integer nuclear spins with degenerate 
ground states of nuclear quadrupoles. 

III. RESONANT SUSCEPTIBILITY 

The main attention of this work is paid to the TLS resonant dielectric susceptibility e res 
which shows the large deviation from the conventional tunneling model for T < 5mK. We 
begin from reviewing the nature of the contribution of tunneling systems to the dielectric 
constant. Although the resonant dielectric susceptibility has the well established behavior 
(for reviews, see [3, 4] and references therein), our analysis will be helpful for understanding 
the effect of the nuclear quadrupole interaction. It is also useful for readers not too familiar 
with the TLS dielectric response. 

Consider a single TLS polarization due to the external electric field F. We suppose that 
tunneling atoms possess a nonzero charge. Then tunneling between the right and left wells 
changes the dipole moment of TLS. The TLS dipole moment operator can thus be expressed 
in the terms of pseudospin (see Eq. (3)) 



Here /x is a dipole moment of a tunneling system. The interaction of the external field F 
with a TLS can be written as 



In general, there can be a contribution to the dipole moment proportional to the off-diagonal 
operator s x . This results from a change of the barrier due to electric field. Employing the 



fl = fis z . 



(16) 



V = -F/xs 2 . 



(17) 
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experimental data [3] and theoretical estimates [24], we can argue that such term leads to 
much smaller effect than that from the "diagonal" term Eq. (16). Therefore we neglect it. 

The external field effect application can be taken into account by introducing the field- 
dependent asymmetry energy 

A(F) = A + F£. (18) 

Thus, the energies E a , a = 1, 2, ...Z, of all Z = 2 • (21+ l) n eigenstates of a tunneling system 
become dependent on the external field. Remember that Z = 2 in the absence of the nuclear 
quadrupole interaction. 

We are interested in a linear response of a TLS ensemble, i.e., response to an infinitesimal 
electric field. Then the dipole moment of the eigenstate a can be expressed as 

dE a dE a 

» a = ~ ~cW = -» ~dK ■ (19) 

The total TLS dipole moment can be expressed as a sum of contributions of all the Z 
eigenstates a, weighed by the Gibbs population factors P a 

z 

K a P. 

<9A,; 



a=l 

The population factor P a is given by the equilibrium distribution for unperturbed TLS's 

P Q (A ,A;T) = _^^f\ N . (21) 



E 7=i exp (--^ 
Finally, the susceptibility of a given TLS is determined by 

X*-W b -^ {-^ bPa ~dA^ + ^dF b ~dK J • (22) 

a=l v ' 

Here the indexes a and b denote the Cartesian coordinates of the corresponding vectors 
and tensors. The first term is associated with the adiabatic excitation of the TLS due to 
the field-induced change in its eigenenergies. This contribution reaches its maximum when 
TLS has an asymmetry A smaller than its tunneling amplitude A . It is called a resonant 
contribution [3]. The remaining term is associated with the change in populations of TLS 
energy levels induced by the external electric field. Such changes take place by the TLS 
transitions between different quantum states. Therefore, this contribution corresponds to 
the relaxational term [3]. 
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For temperatures T < 50mK and an external alternating field F with frequency v > 
100Hz, the relaxation contribution is negligibly small [3, 5, 15] because a relaxation rate 
of TLS populations becomes much smaller than the field oscillation rate v. Therefore, the 
relevant TLS's cannot adjust their thermal populations to the rapidly changing field. In 
the range of interest T ~ 5mK we can neglect the relaxational contribution and restrict our 
consideration to the first term on the right-hand side of Eq. (22). 

The contribution of a single TLS should be summed over all TLS's belonging to the 
system. This is equivalent to averaging the susceptibility Xab in Eq. (22) over energies, 
tunneling amplitudes, dipole moments, and other possible relevant parameters. Averaging 
over the directions and absolute values of TLS dipole moments is straightforward and we 
can rewrite the resonant TLS contribution in the form 

X = ^£^(Ao,A;T)x a (A,A )^, 

Xa(A,Ao) = -|^. (23) 

Here /xjj is the average square of the TLS dipole moment. Also, we have ignored the corre- 
lations between TLS dipole moments with the other tunneling parameters. This is justified 
by the available experimental data [3]. The average < ... > implies the integration of the 
single TLS response over dAdA / A in accordance with the postulated TLS distribution 
Eq. (2). 

Let us first consider the behavior of the dielectric constant in the zero-temperature limit. 
In this case it suffices to take into account the ground state only in Eq. (22), calculating 
the resonant dielectric constant 



The upper integration limit W represents the maximum tunneling amplitude, while the 
integration limits for an asymmetry parameter A are set to ±oo because its absolute value 
can be much larger than the tunneling splitting A . The integral over asymmetries can be 
evaluated as 

6E g (A ,A) *=°° 



dA 
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For a large asymmetry |A| 3> A , the ground state is determined by the minimum energy 
state of a particle in the left A > (see Fig. 1) or the right well A < 0. So, for large A the 
ground-state energy behaves as E g — | A | /2, while the tunneling amplitude and nuclear 
interaction give rise only to the small correction. Then the expression in Eq. (25) equals 
unity. Therefore, Eq. (24) results in the divergent integral 

3 i„ A„ 1 ' 

It is remarkable that this result does not depend on whether or not the quadrupole interaction 
exists. So, the question remains: what happens at finite temperature? 

The answer is clear for TLS's without a nuclear quadrupole interaction. In this case each 
TLS has only the two states, i.e., the ground state g and the excited state e with energies 
E g>e = =f(1/2)a/A 2 + A% (Eq. (1), cf. Ref. [3]), respectively. The susceptibilities of excited 
and ground states differ by the sign since 

d 2 E e = = 1 A 2 

OA 2 OA 2 2(A2 + A2) 3/2 ' 

The response comes mainly from resonant TLS's having | A| < A . Therefore, both levels of a 
TLS with A < T are approximately equally populated. For this reason, the contribution of 
the excited state to the dielectric constant cannot be neglected. Moreover, the contributions 
from the excited and ground state nearly cancel each other if A < T. 

Calculating the finite temperature resonant susceptibility, one should consider only TLS 
ground states and cut-off the integral Eq. (24) at the lower limit given by A ~ T. This 
leads to the well-known logarithmic temperature dependence 

This result is valid as long as the temperature exceeds the energy of the quadrupole 
interaction nb. Next, we discuss the case T ^ b. TLS's with small tunneling amplitudes 
A < nb still contribute to the resonant susceptibility. They can be represented by pairs 
of lowest nuclear quadrupole levels in the right and left wells because the higher levels are 
separated by the gap b 3> T from these two lowest ones. They are coupled with each other 
by the tunneling amplitude A reduced by the overlap factor Eq. (11) (see Fig. 2), i.e., 



A oV n . (29) 
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FIG. 2: Two level configuration with splitted energy levels. 



These two lowest levels can be treated as a new TLS. Since only TLS's with Ao* > T 
contribute to the permittivity (see Eq. 28), this defines the renormalized lower cut-off 



^01 



Trf n . (30) 



Substituting this cut-off into integral Eq. (28) yields in the limit T — > 

* - p 4 C t £ = ^ iMW/T) - - ,n(1/ " )) • (31) 

Thus, due to the quadrupole interaction, we predict qualitatively a noticeable reduction of 
the TLS contribution to the dielectric constant at low temperatures. We believe that this 
reduction can explain the plateau in the temperature dependence of the dielectric constant. 
It follows from the above analysis that in the whole energy interval 

A ,J A »' A ° >>n6 ' (32) 
[ A n n , A < nb. 

This renormalization results in a gap in the distribution of the effective tunneling amplitude 
A * of tunneling systems 

0, nbr] n < A * < nb (33) 



P(A *) = < 



p 

A , ' 



A * < nbr] n . 



Using the latter result to estimate the dielectric constant temperature dependence in the 
expression similar to Eq. (31), one can obtain the plateau in the temperature dependence 
of the dielectric constant within the range nbrf 1 < T < nb. At T > nb one should use 
the standard tunneling model result Eq. (28), while at T < nbrf 1 the resonant dielectric 
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constant obeys Eq. (31). Thus, our qualitative arguments can explain the behavior observed 
experimentally. 

In the following part of this paper we investigate the two regimes of Eqs. (28), (31) and 
the crossover between them with the higher accuracy. 

IV. PERTURBATION THEORY APPROACH 

The renormalization of a tunneling amplitude A affects the thermal and kinetic prop- 
erties of the TLS ensemble and provides a minimum energy splitting of TLS's having zero 
asymmetry A = 0. In particular, the resonant susceptibility is determined by similar res- 
onant TLS's with small asymmetry |A| < A . Therefore, we will study the most relevant 
case, i.e., A = 0. The resonant permittivity we are interested in here is determined by the 
ground state of the tunneling systems with the energy difference between the ground and 
first excited states larger than the thermal energy. In this case the contribution of higher 
excited states is insignificant because due to an exponentially small probability of their oc- 
cupation. A positive contribution of a ground state to the dielectric constant reflects the 
general fact that the ground state minimizes the energy of the system. Thus, the external 
electrical field aligns the ground state dipole moment along the field direction. The suscep- 
tibility of excited states can be negative like in the case of Eq. (27) and also in the case of 
TLS's is affected by the nuclear quadrupole interaction as we will show below. Therefore, 
when the temperature becomes comparable or larger than energy splitting between ground 
and excited states, TLS susceptibility becomes negligible. Thus, a structure of the ground 
state is only important for the resonant susceptibility of the system. 

Next, we consider the effect of the quadrupole splitting on the ground state of the tun- 
neling system. In the case of vanishing asymmetry energy A = the Hamiltonian of the 
tunneling system can be represented in the form 



The total wave function of a tunneling particle is a product of the coordinate (pseudospin) 
wave function and the nuclear one. 

If the tunneling amplitude A is large, the ground state is described by the wave function 
\ghyb) whose coordinate part is symmetric and the tunneling particles are shared equally 




(34) 
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between the two wells. In this case one has < s x 1/2, < s z and the energy of the 
nuclear quadrupoles is given by the "mean" Hamiltonian 



H r + Hi 



(35) 



This regime is called the hybridized one. 

In the opposite limit of the small tunneling amplitude, i.e., in the localized regime, one 
can neglect the tunneling term s x in Eq. (34) and the energy minimum corresponds to either 
the ground state of the Hamiltonian H[ or to the the ground state of the Hamiltonian H r . 
Then, the influence of the tunneling s x term on the ground state is insignificant due to the 
strong reduction in the effective tunneling amplitude (see Eq. (32)) because of the small 
factor r] n . Let us find the crossover between the hybridized and localized regimes. 

In the hybridized regime the ground state energy is 



A 1 

— + - \ Qhyb 



H r + Hi 



9hyb 



(36) 



while in the localized regime it is 



9r 



H r 



9r 



9i 



H 



9i 



(37) 



Let us introduce a parameter A* describing the ground-state quadrupole energy difference 
between the two limiting regimes. It is the reorganization energy corresponding to the 
transition from the hybridized state to the localized one 



A* = - (gi 



H 



9l) + ( Qhyb 



H r + Hi 



9hyb 



(38) 



Comparing Eqs. (36) and (37), one finds that the hybridization regime is realized if 



A > 2A„ 



while in the localized regime this inequality changes its sign, i.e., 



A < 2A*. 



(39) 



(40) 



Let n be the number of atoms of a TLS experiencing nuclear quadrupole interaction. The 
parameter 

6* = A*/n (41) 
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represents the reorganization energy per a tunneling atom. Let us estimate the parameter 
6* for the case when the quadrupole interaction is described by Eq. (5). One has 

£i = *(/ x 2 -/(/+l)/3), (42) 
H r = 6 ((4 cos + sin 0) 2 - /(/ + l)/3) . (43) 

Then one can calculate the parameter 6* by using Eqs. (38), (42), (43). To be more specific, 
consider the case b > and integer spin I = 1,2. The results of the numerical calculation 
of a single atom parameter in the model Eqs. (42), (43) for different rotation (mismatch) 
angles are represented in Fig. 3. One can see that the reorganization energy is always 
comparable with the nuclear quadrupole interaction energy. 




12 3 



Rotation angle (j)(rad) 

FIG. 3: Effective interaction energy 6* of nuclear spins vs. angle between the nuclear quantization 
axis in two wells. 

In the hybridization regime Eq. (39) the quadrupole interaction is a weak perturbation. 
Then, the multiplet structure proves to be insignificant when one calculates the contributions 
of these tunneling systems to the susceptibility. In other words, the latter can be calculated 
within the standard TLS approach. 

However, in the localized regime Eq. (40) a multiplet structure is the decisive feature. 
In this regime a renormalization of the tunneling amplitude (32) becomes important. If a 
tunneling amplitude A is small enough, the ground state will be a superposition of the 
ground states in the isolated left and right wells (see Fig. 2). Let us estimate the upper 
limit of A below which the two-level approximation for the ground state is still valid. 

The nuclear-spin ground state in each well can be treated as non-degenerate (see. Sec. 
II B). One can approximately construct the ground state and the lowest excited state of the 
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tunneling system as a superposition of the unperturbed ground states in the two wells. The 
contrubution of higher excited states are neglected. This is justified by calculating their 
contribution to a lowest-order perturbation theory. We estimate the correction factor c for 
the ground state amplitude as 



where i labels all states of the right well. 

We determine the parameters of the regime for which the second term can be neglected as 
follows. The lowest excited states in each well are separated from the ground state by some 
characteristic energy b ~ fkj where u is a frequency of the nuclear quadrupole resonance. 
The next group of states is separated by the energy gap ~ 2b. The vast majority of states 
have energies exceeding that for the ground state by the energy A ~ nb. Because of the 
large statistical weight, these states provide the main contribution to Eq. (44). Therefore 
we may replace the denominator in Eq. (44) by A. Then, the sum of the overlap factors 
|< gl\i >| 2 in the numerator can be well approximated by unity. The requirement that the 
second term in Eq. (44) is small results in the condition 



This condition is weaker than Eq. (40). Therefore, c ~ 1 is justified when (40) is fulfilled 
and one can ignore the excited states in both wells. 

Thus we can specify two domains of system parameters. One is given by large tunneling 
amplitudes A > 2A*, where the nuclear quadrupole interaction can be ignored. The other 
is given by small tunneling amplitudes A < 2A*. In the latter case we can restrict our 
consideration to the pair formed by the two lowest energy states. These eigenstates are 
superpositions of the ground states in the left and right potential wells. A non-perturbative 
approach developed in the next section, where nuclear spins are replaced with oscillators, 
leads to the similar results only with some minor deviations. 

At this point we turn to the analysis of the dielectric constant using the structure of 
TLS ground states described above. For an improved analysis, it is convenient to write the 




(44) 



A < 2A. 



(45) 
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resonant dielectric constant Eq. (22) in the following form 



X = y E ^ ( A °> A - T ) Xa (A = 0, A)) 



+^ E (( P « ( A °> A ' T ) " P « ( A °> A - °)) x (X- ( A o> A) - Xa (A = 0, A))) 



3 



a 



+^ E ( A °> A - °) ( A o> A) - Xa (A = 0, A))) . (46) 

a 

The logarithmic temperature dependence, often serving as an evidence for the TLS effects [3], 
comes entirely from the first term. This is is due to logarithmically uniform distribution of 
the tunneling amplitudes Eq. (2). The logarithmic divergence of the first term is suppressed 
for TLS's with the energies of the order of the thermal energy due to the compensation of 
a positive contribution from the ground state by negative contributions of excited states 
(see below). The third contribution is a temperature-independent constant. We will ignore 
it as a background correction to the susceptibility which does not affect its temperature 
dependence. 

In the absence of the nuclear quadrupole interaction the important first and second terms 
can be evaluated. One can reexpress the first "logarithmic" term as 

Xw = ^/ W ^tanh(A /(2T)), (47) 

and the second "thermal" contribution as 



rdAof 
Jo A J 



Xt =^21 / _l dA\ tanh 



2T 



A ° s/2 - 5(A) I . (48) 
^2(A 2 + A2) 3/2 

In the second integral the upper cut-off for the tunneling amplitude is replaced by oo because 
the integrand decreases exponentially at large energies. 
The first contribution can be written as 



Xu » = -3-\ ln \2f)- 11 
f°° \n(x)dx 

h ~ / -^TTTT ~ "0.82. (49) 
Jo cosh {xj 
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The "thermal" contribution is given by the dimensionless integral 

-+oo 




/ °° °°^( tanh (v / ^+^- 1 )) 

- <%) ) . (50) 



x 2 



2 ( X 2 + y 2)3/2 



It vanishes because the upper limit of the integration over x has been replaced by 00. This 
can be demonstrated by using, for instance, the trigonometric substitution x = a cos(/3), y = 
asm((3) and evaluating the integral over (3 first [25]. 

Below we calculate the susceptibility of tunneling systems for different temperature 
regimes. According to the above analysis (see Eqs. (39), (40)) it is convenient to divide all 
tunneling systems into two parts depending on either Ao > 2 A* or Ao < 2A*. 

In the first case A > rib* the effect of the nuclear quadrupole interaction is negligible 
and one can use the standard two-level approximation in order to estimate the contribution 
from these TLS's into the susceptibility as 

/^f ^tanh^%% (51) 



J2K A V 2T 7 

Consider temperatures T > 2A*. With E(A Q ,0) = A one can estimate this integral as 
follows (cf. Eq. (49)) 

W A* 

/i« In— --^ + 0.82. (52) 

The second term on the right-hand side is small in comparison with the first one and can 
be omitted. If the temperature is low T < 2A*, the tangent in Eq. (51) equals unity and we 
have 

W 

h = hx—. (53) 
nb* 

Now consider the second group of tunneling system for which Ao < A*. For them the 
multiplet structure becomes important and we can confine ourselves to ground state levels 
of the multiplet in each well. These levels are coupled by the tunneling amplitude A ^ n . 
Then one has 

Assume that T < 2A^ n . Then 

h ~ In — . * . (55) 
T/r] n K J 
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This expression should be added to the term Ii given by Eq. (53), also giving a contribution 
the temperature region considered. Thus, the total integral is given by 

Ji + / 2 = ln^U (56) 

In the intermediate temperature range 2A* > T > 2\^r] n Eq. (54) is inapplicable strictly 
speaking. However, we expect that the contribution to the dielectric permittivity within 
the range 2A* > A > 2\*r] n is small because this range corresponds to the gap in the 
distribution of effective tunneling amplitudes (Eq. (33)). In addition, the lower limit of the 
low temperature dielectric constant Eq. (56) coincides with the upper one of the dielectric 
constant in the high temperature range Eq. (52). Thus, in the whole temperature region 
one has 

e res w (P/i 2 /3) HW/T), 2K < T, 

e res « (P/i 2 /3) hi(W72A*), 2A*??" < T < 2A*, 

e res w (P/i 2 /3) ln(WWT), T < 2A^ n . (57) 

Our perturbation theory analysis is a phenomenological one and cannot be exact because 
we do not know the nuclear spin Hamiltonian. However, it is expected that the main qual- 
itative features of the systems behavior are reproduced. For a more quantitative and non- 
perturbative analysis given in the next section, we replace the nuclear spins with oscillators 
and investigate the tunneling amplitude behavior in such a toy model. 

V. TOY MODEL 

Since the nuclear spin interaction in the right and left wells is unknown, we consider in- 
stead a toy model in which nuclear spins are replaced by classical oscillators. This "bosoniza- 
tion" approach to nuclear spins is justified when one deals with the low energy states of many 
spins. In particular, it has been used by Prokofev and Stamp to investigate the nuclear spin 
interaction effect on the large spin tunneling [26]. We consider a symmetric TLS (A = 0) 
characterized by a coherent tunneling amplitude A and coupled to n oscillators representing 
the nuclear spins of n tunneling atoms forming the TLS concerned. All oscillators have the 
frequency f2 and the mass M. These n oscillators are linearly coupled to the TLS and have 
the shifted equilibrium positions X; L = ±x /2, i — 1, ..n when the TLS occupies the right or 
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left wells, respectively. Then the Hamiltonian of the system can be expressed as 

i=i v 

n 

- s z MQ 2 x ^2xi. (58) 
i=i 

The spin values s z = ±1/2 stand for the TLS residing in the right or left wells, respectively. 
The models similar to Eq. (58) have extensively been studied within the polaron theory 
(See e. g. Refs. [27, 28] and refrences therein). Following the standard approach one can 
determine an approximate ground state of the problem by minimizing the "classical" part 
of Eq. (58). This excludes the kinetic energy term of the oscillators. The "classical" energy 
of the system can be expressed as 



E cl - -\ 



Ag+(Mn'soS>) (59) 



where the spin-Hamiltonian has been replaced by its ground state energy, thereby using 
the relationship E ground = — |a/A 2 + Aq for the Hamiltonian — As z — A s x . Oscillators are 
treated classically. This is justified when their number is large, i.e., n >> 1. We suppose 
that this regime is applicable here. 

The total energy can be minimized with respect to the oscillator coordinates Xj. Then, 
one has for the derivatives of the classical energy Eq. (59) with respect to all coordinates 

dE d (Mn 2 x ) 2 X 
= — - = / ° J = + MQ Xi, 

dXi 2^/A 2 + (Mn 2 x X) 2 

n 

X = J2^- (60) 
i=i 

These equations can be solved analytically by taking into account that at the energy mini- 
mum all equilibrium coordinates Xi are identical, that is X{ = X/n. The following relations 
are found 

A > 2A* : X = 0, 



A* = nh, K = MQ 2 x 2 /2. (61) 
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The quantity K has been introduced to describe the change in the oscillator energy induced 
by its interaction with TLS. For spins, this energy is equivalent to the quadrupole splitting 
energy 6* (see Eq. (41)). 

Thus, depending on the relationship between the oscillator energy A* and the tunneling 
amplitude A , a TLS ground state can have different structures. When tunneling is stronger 
than the TLS interaction with oscillators, i.e., when A > 2A*, the tunneling atoms are 
distributed equally between the two wells and all oscillators have their energy minimum at 
Xi — 0. This state is energetically most favorable because an identical internal structure for 
the left and right states minimizes the energy of the system. In other words, the tunneling of 
TLS is so fast that the oscillators see it in its average position in both wells simultaneously. 
This state is equivalent to the hybridized state considered previously in Sec. IV. 

In the opposite case of weak tunneling A < 2A* the symmetry between the right and 
left wells is broken because of a strong displacement of oscillators. In this case the system 
has two energy minima depending on the sign in the definition of the displacement X in 
Eq. (61). For the nuclear quadrupole interactions, these two states are represented by the 
nuclear spin configurations minimizing the energy of TLS's localized either in the right or 
left wells. 

Both the energy minima are still coupled by tunneling, but the tunneling amplitude 
is much weaker then in the case of vanishing oscillator displacements. Consider the non- 
adiabatic tunneling regime applicable to the nuclear spins [26] . In this case we may express 
the effective tunneling amplitude A * as the product of the coherent tunneling amplitude 
A and the overlap integral < l\r > of the left and right states of the environment, i.e., 

A * = A -(Z|r). (62) 

In order to estimate the overlap integral, one can use the harmonic approach. We will use it 
for the wavefunctions in the left and right wells and consider the case of zero temperature, 
reasonable if T < 2A*. The domain T 2A*, where excited states are important, is relatively 
narrow at the large number of oscillators n and can be approximately ignored due to a 
logarithmically weak dependence of the dielectric constant on T. 

In the harmonic approach one may expand the energy Eq.(59) near the local minima 
given by Eq. (61) up to the second order in a coordinate displacement. This expansion can 
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be written as 



2 

hi 



d 2 E, 



4,=^=-si«,-^i i. 

where 5ij is the Kronecker symbol. The expansion is written near the potential minimum at 
the right well, while for the left well the sign of X should be changed to the opposite. 

The harmonic part of the Hamiltonian Eq. (63) can be represented by n independent 
modes including the symmetric mode 




i=i 

and n — 1 asymmetric, degenerate modes a 



« = ^f>, n s = n Ji-^, (64) 



a = Y,$Xi, V a = n , = ^|c"| 2 = l. (65) 



i=l i i 

A factor 1 / y/n is introduced in Eq. (64) to conserve the commutation rules between the 
coordinate and momentum operators. Only the symmetric mode interacts with the tunneling 
motion between the right and left minima, while asymmetric modes remain unchanged and 
can therefore be ignored. Therefore the overlap integral in Eq. (62) is given by the overlap 
between two ground-state wavefunctions of the symmetric harmonic modes with the mass 
M and the frequency Q s , and the equilibrium positions shifted by ±\X\/y/n Eq. (61) from 
the origin. Using the Gaussian wavefunctions for oscillator ground states, one can express 
this integral in the form 

. MttX 2 ' 

{I I r) — exp 

=« »{-^[ 1 -{t;)) )■ (66 » 

This result can be used to characterize approximately the tunneling of a TLS consisting of 
n atoms coupled to the nuclear spins. The overlap integral for a single oscillator % is 

v = (U | n) = exp ( ^ ) (67) 
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in analogy to the single-spin overlap integral introduced before. The parameter A* can be 
represented by the single atom quadrupole splitting energy multiplied by the number of 
atoms in a TLS, i.e., by 

A* ~ nb*. (68) 

Thus we may use the following approximate relationship between the initial (coherent) 
tunneling amplitude A and the effective tunneling amplitude A * 

3/2 N 



*0* 



A exp f-n 111(1/77) ^ - (^) ) J " ( 69 ) 



In the following we will use this relationship in order to describe the dielectric response of 
TLS. 

The tunneling amplitude distribution modified by Eq. (69) has a dip (pseudogap) due 
to sharp decrease in A * when the coupling strength A becomes smaller than the effective 
rearrangement energy A* = nb*. We can define the modified distribution of TLS's over Aq* 
for A * < A* by using 

P(Ac) oc ^ x 5 |ao* - A exp (-n ln( V ) - (^j j j j (70) 

(TD 



A °*l + 3nln(l/,) (%)yi-(%) 2 



where the amplitude Ao is an implicit function of Ao* defined by Eq. (69). In the case of 
a large tunneling amplitude Ao* > A* the distribution is logarithmically uniform because 
A = A * in that regime. 

For the proof of Eq. (71) we employ 

5(a-<f>(x)) = ^—5(x-x), (72) 

where x is the solution of equation a — (x) = O.We identify 

a = A *, x = A (73) 

x N 2 \ 3/2N 

(x) = x exp j — n \n(rj) [ 1 — ( — ) ) | (74) 

A* 
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and derive 



/ (A o )=exp f-nlnfa) ^1-0^ 



3/2 N 



x(l + 3nln (l) )|l-[^) 2 )' /2 (x 0S2 



(75) 



By substituting Eq. (69) into (75), we obtain 



'/ (A ) = ^ | l + 3nlnfa) ( 1 



A \^ Mo w 



A* 



(76) 



When this expression is put into Eq. (72) and then into (70), we obtain 



dA 1 

^ l + 3nln(,)(l-(t) j (!f 



5 A - A 



A =A 



which coincides with Eq. (71). The modified distribution of tunneling amplitudes 
W(ln(A „) = P(A *)A * is shown in Figs. 4, 5. 
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FIG. 4: The dip in the logarithmic TLS distribution over tunneling amplitudes due to the nuclear 
quadrupole interaction. The logarithmic density of states P(Aq)Aq is shown for / = 1 and different 
numbers n of tunneling atoms. The case n = corresponds to the lack of the interaction. The 
rotation angle is set to 1.23 rad. 



We have chosen a representative value of = 1.23 rad for the rotation angle (see Eq. 
(9)) in a field of axial symmetry bl 2 a and b > 0. In accordance with Eq. (9) and Fig. 3 the 
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FIG. 5: The dip in the logarithmic TLS distribution over tunneling amplitude due to the nuclear 
quadrupole interactions. The logarithmic density of states P(Ao*)Ao* is shown for I = 2 and 
different numbers n of tunneling atoms per TLS. 

effective interaction parameters are r\ = 1/3 and b* = 0.33mK for 1=1, and i] = 1/3 and 
5* = 0.85mK for 1 = 2. The nuclear quadrupole effect shows up at an energy which is about 
three times larger for / = 2 then for I = 1. 

One can see from Figs. 4, 5 that the dip in the distribution of tunneling amplitudes 
appears in the domain (A* exp(— nln(l/r/)), A*). This becomes sizable for n = 4. A similar 
behavior is found for the other values of the overlap integral rj. This dip causes a change of 
the temperature dependence of the resonant susceptibility as discussed below. 

The resonant dielectric susceptibility of the TLS ensemble can be estimated from the 
"logarithmic" integral 



By integrating the distribution with the dip we obtain the results shown in Figs. 6, 7. 

The quadrupole interaction parameter b = lmK and the rotation angle ~ 1.23rad 
are the same as before. One notices a plateau in the dielectric constant separating two 
logarithmic temperature dependences characterized by the same slope at sufficiently large 
number n of atoms per single TLS. In the experiment [15, 16] a plateau is found for T < 5mK, 
while the low temperature edge of the plateau has not been observed because it requires 
much lower temperature than those reached experimentally. A temperature of 5mK gives a 
reasonable estimate of the quadrupole interaction energy. 




(77) 
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FIG. 6: Resonant dielectric constant affected by the nuclear quadrupole interaction in the case of 
the nuclear spin 1 = 1. 




Temperature, T (mK) 



FIG. 7: Resonant dielectric constant affected by the nuclear quadrupole interaction in the case of 
the nuclear spin 1 = 2. 

In all cases the plateau begins at the characteristic maximum temperature 

and extends down to the characteristic minimum temperature 

T mm ~ \* V n . (79) 

For T < T m in the behavior of the standard tunneling model is restored and the dielectric 
constant shows a logarithmic temperature dependence with the same slope as at high tem- 
peratures T > T max . The logarithmic width of the plateau In (T max /T min ) is thus directly 
proportional to the number of atoms per single TLS. 
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We expect that in the plateau regime the dielectric constant can be controlled by the 
external magnetic field. In fact, the overlap integral of the nuclear spin states in the left and 
right well can be affected by the magnetic field. This effect is discussed in the next section. 

VI. EFFECT OF THE EXTERNAL MAGNETIC FIELD ON THE ANOMALOUS 
DIELECTRIC PROPERTIES AT ULTRA-LOW TEMPERATURES 

The application of the external magnetic field affects the orientation of nuclear magnetic 
moments. When Zeeman splitting becomes comparable with the nuclear quadrupole inter- 
action, the quantization axes in both potential wells of the given TLS are aligned with the 
direction of the magnetic field. Accordingly, the mismatch of nuclear quadrupole axes will 
be reduced by the field and the overlap integral of nuclear quadrupole states in different 
wells increases and approaches unity in the high field limit. The effective interaction con- 
stant b* vanishes in that case. Thus, high magnetic field reduces the influence of the nuclear 
quadrupole interaction and the dielectric constant behaves like in the standard tunneling 
model. In this section we estimate the effect of the external field on the low-temperature 
resonant dielectric constant by using the simple oscillator and axially symmetric models for 
the nuclear quadrupole interaction formulated in previous sections. 

Consider the effect of the external magnetic field on parameters characterizing the mis- 
match of the nuclear spins in the right and left wells. The set of parameters necessary to 
characterize the given TLS includes the overlap integral for single nuclear spin 77, the num- 
ber n of nuclear spins involved into the tunneling process and the characteristic nuclear spin 
interaction energy b* (see Eq. (69)). The number n is field-independent, while the overlap 
integral rj and the single atom interaction constant 6* are subjected to the changes with 
the external field. We investigated numerically the change of the parameters 77 and b* in 
an external field for the axially symmetric quadrupole interaction. The nuclear spins are 
1=1 and 1 = 2 with the same angle cos(0) ~ 1/3 and the nuclear quadrupole interac- 
tion constant b = lmK as has been assumed in previous sections. The sign of the nuclear 
quadrupole interaction constant b > is chosen as previously. Therefore, the ground state 
has the zero spin projection onto the quantization axis. The zero-field dielectric constant for 
these regimes is shown in Figs. 6, 7. We have calculated numerically the dependence on the 
external magnetic field B of the overlap integral between the nuclear spin ground states in 
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FIG. 8: The effect of the magnetic field on the overlap integral. The lmK Zeeman splitting 
approximately correspond to field 5T. 
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FIG. 9: The effect of magnetic field on the effective single-atom interaction constant 6* 



the left and right wells. Also we have estimated numerically the dependence of the effective 
interaction constant on the magnetic field. The calculations have been made as follows. 
The quadrupole Hamiltonian is chosen in the form b{I a ) 2 with a the left-well axis a = x 
and the right-well axis rotated by the angle = 1.23rad in the x — y plane with respect to 
the x— direction. Then the Zeeman term describing the interaction of the nuclear magnetic 
moment with magnetic field, is augmented to the Hamiltonian (5). The magnetic field with 
a fixed absolute value B has been generated in the random directions. The overlap integral 
of ground states r\ as well as the effective interaction constant 6* have been computed and 
then averaged over ~ 10 4 realizations of the random field. The results of calculations are 
shown in Figs. 8, 9. 

We express the magnetic field in units of the Zeeman energy splitting (mK). In fact, it 
is unclear which nuclear spins are most important in the experiment. A reasonable scale we 
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FIG. 10: Effect of a magnetic field on the dielectric constant in the case of spin 1 = 1 and the 
number of atoms per TLS n = 4. 
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FIG. 11: The effect of magnetic field on the dielectric constant in the case of spin 1=1 and the 
number of atoms per TLS n = 4. 



use for estimates is about lmK spin splitting in a field of B ~ 5T. 

As is clear from Figs. 8, 9, the overlap integral increases monotonously with the field 
and approaches unity when \iB ^> lmK. The effective interaction constant 6* vanishes in 
the same limit. Therefore, the application of magnetic field results in the disappearance 
of the plateau in the temperature dependence of the dielectric constant and the standard 
logarithmic dependence is restored. To examine this effect, we compute the temperature 
dependence of the TLS resonant dielectric constant at various fields. We have used Eq. (77) 
with Eq. (69) for the effective tunneling amplitude and input parameters rj and 6* obtained 
from our field-dependent calculations (Figs. 8, 9). 

The field-dependent dielectric constant is shown in Figs. 10, 11 for n — 4 atoms per single 
TLS and nuclear spins 1 = 1 and 1 = 2, respectively. The field effect is similar in both cases. 
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However, a stronger field is needed to suppress the nuclear quadrupole interaction in the 
case of larger spin 1 = 2. 

We can interpret the field dependence as follows. A relatively small field [iB < 0.2mK 
and B < IT affects the resonant dielectric constant only for low temperatures T < 5mK. 
At higher temperatures the field effect can not be viewed beeing a small correction. The 
further increase of the field leads to linear reduction in the width of the plateau from both 
the high- and low-temperature sides. When the field is sufficiently high fxB > lmK, and 
B > 5T, it competes with the nuclear quadrupole splitting and the plateau narrows by 
orders of magnitude and almost vanishes at fiB ~ 4mK. Then the standard tunneling model 
behavior is restored. It would be very interesting to perform measurements of the dielectric 
constant at very low temperatures, i.e., for T < 5mK and a strong magnetic field around 
10T. In this case the plateau in the dielectric constant should disappear. We hope that our 
work will stimulate such measurements. 

VII. RELATIONSHIP OF THE MODEL TO REAL SYSTEMS 

In this section we address two questions. First question is whether our choice of param- 
eters for the nuclear quadrupole interaction is justified and these parameters can be used 
to fit the existing experimental data. The second question is related to glasses having no 
(or negligable) nuclear quadrupole interaction. These glasses (mylar, SiC^) should not show 
any anomalies in the dielectric constant temperature dependence. The experimental data 
related to these glasses are discussed below. 

According to the previous consideration we expect that the plateau observed experimen- 
tally in the temperature dependence of the dielectric constant for T < 5mK is due to the 
fact that the parameter nb (see Eq. (6)) amounts to a value of 5mK. Experimentally, a 
qualitative change in the TLS resonant dielectric susceptibility is seen in the temperature 
range of T ~ 5 — lOmK (see Refs. [13, 15, 16]). There is no problem to intepret the satura- 
tion behavior using our formalism within the experimental accuracy choosing properly the 
fitting parameters [21]. However the question arises whether or not the nuclear quadrupole 
interaction is sufficiently strong to display itself in this temperature range. We also expect 
that glasses lacking the nuclear quadrupole interaction should not show any deviations from 
the standard model. 
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The nuclear quadrupole interaction constant b is defined by Eq. (6). The values of the 
electric field gradient are close to each other in different glasses, while the nuclear quadrupole 
moments Q differ strongly for different elements, leading to the broad distribution in the 
observed energies of the nuclear quadrupole resonance hwo ~ b Ref. [29]. The data for the 
chemical elements present in glasses for which the low temperature measurements have been 



made are summarized in Tab 


e 1. 


Nucleus 
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Hn 


Q(barn) 






1/2 


4.837 








D= 2 H 


1 


1.213 


0.00286 


0.7 • 10~ 2 Ref. [18] 


12 C 














16Q 














23 Na 


3/2 


2.863 


0.104 


0.14 (Na), 0.4 (NaF) Ref. [30] 


27 A1 


5/2 


4.309 


0.147 


0.9 in Al Ref. [31] 


29 Si 


1/2 


-0.962 








29 K 


3/2 


0.505 


0.0585 




135 Ba 


3/2 


1.82 


0.160 


0.86 (BaBi0 3 ) Ref. [32] 



Table 1. Nuclear spin, magnetic moment, quadrupole moment and frequency of nuclear 
quadrupole resonance for different chemical elements, possibly participating 
in the tunneling. 

As one can see from Table 1, the typical range of nuclear quadrupole interactions for Na, 
K, Al, Ba nuclei is around but slightly less than lmK. This interaction should be larger 
in glassy materials because the above-mentioned elements are bound with non-metal atoms 
there. The covalent and ionic bonds are expected to be stronger than metallic bonds, in 
particular, due to the absence of screening. Therefore, their binding energy and, accordingly, 
the electric field gradient affecting the nuclear quadrupole interaction are expected to be 
larger in these glasses. 
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Glass 


Nuclei 


Xl„/ fmrO 

sat v / 


Exp. 


Refs. 


Mylar fCinH s 4 )„ 


No 


< lmK 


£ 


Ref. [15] 


5%K-Si09 


K 


4mK 


£ 


Ref. [15] 


10%K-SiO9 


K 


4mK 


£ 


Ref [151 


BK7 [33] 


Na 


5mK 


£ 


Ref. [15] 


SiO* 


No 


8mK 


£ 


Ref. [15] 


BaO-Al 2 3 -Si0 2 


Al,Ba 


~ 5mK 


e 


Ref. [16] 


a— Si02 


No 


< 2mK 


V 


Ref. [34, 35] 


BK7 


Na 


< 5mK 


V 


Ref. [36] 



Table 2. Saturation temperature below which the dielectric constant e and/or sound 
velocity v become temperature-independent. 

The experimental data indicating the strong changes in low temperature dielectric and, 
probably, acoustic properties are summarized in Table 2. The saturation in a temperature 
dependence of the dielectric constant below the temperature T sat takes place in all materials 
containing Na, K, Al or Ba which have relatively high quadrupole moments (see Table 1). 
The high saturation temperature observed in BaO-Al 2 3 -Si0 2 is most likely due to large 
values of the Ba and Al quadrupole moments compared with the other elements in Table 1. 

The absence of the low temperature saturation in the dielectric constant of mylar is in the 
full agreement with the theory. Indeed, mylar is an organic polymer composed of C, H and 
O atoms, for which the most stable isotopes have vanishing nuclear quadrupole moments. 
Similarly according to the most recent experimental data [34, 35] there is no saturation in 
the temperature ependence of the sound velocity in a— Si0 2 having no nuclear quadrupole 
interaction. A saturation in the low temperature dependence of the dielectric constant has 
also been observed in SiO x which poses a puzzle. One possible explanation is that unpaired 
electrons may be present in this material [37]. They act like nuclear quadrupole moments 
in that aspect. 

The problem of interpretation of the acoustic experiment [36] is a knotty one. If the same 
two-level systems control the acoustic and dielectric behaviors of the glasses, the saturation 
should take place below the same T sat both for the dielectric constant and for the sound 
velocity [3]. On one hand, there is no saturation in the logarithmic temperature dependence 
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in a— S1O2 having no nuclear quadrupoles in agreement with the theoretical expectations. It 
is difficult, however, to understand the absence of any deviations from the standard tunneling 
model for the sound velocity measurements in BK7 down to 5mK Ref. [36]. Although the 
major contribution to the dielectric constant and the sound velocity can be due to different 
subsets of TLS's possessing either larger dipole moments or stronger elastic coupling to 
lattice vibrations [38], there must be some amount of TLS's containing quadrupole nuclei 
and contributing to both dielectric and acoustic properties. These TLS's should lead to an 
anomalous acoustic behavior. One possible explanation of the absence of any saturation 
effect is that the measurements have been made for T > 5mK, while T sat 5mK (see Table 
2). An extension of acoustic measurements in BK7 to lower temperatures should help to 
clarify the puzzling situation. 



VIII. DISCUSSION AND CONCLUSIONS 



In this work we have considered various aspects of the effect of the nuclear quadrupole 
interaction on the low temperature properties of glasses. The significance of this interaction 
has been pointed out by Wiirger, Fleischmann and Enss [17]. In the present paper this inter- 
action has been employed to characterize the resonant dielectric susceptibility of amorphous 
solids at ultra- low temperatures T < 5mK where major deviations from the predictions 
of the tunneling model are observed. The analysis of the typical interaction parameters 
shows that the nuclear quadrupole interaction is strong enough to explain these deviations. 
To our knowledge the materials having no nuclear quadrupole interaction and no unpaired 
electrons, i. e. the organic polymer mylar [5], and a-Si02 Ref. [34, 35] do not show the 
deviations from the logarithmic temperature dependence of a dielectric constant or a sound 
velocity down to the lowest temperatures accessible experimentally around T ~ lmK. These 
results support our theory. 

Our theory uses the fact that the nuclear quadrupole interaction is generally different in 
the right and left wells of a two-well tunneling system. Therefore, it affects the coherent 
coupling between the wells. The overall quadrupole interaction effect is governed by the 
relative magnitude of two parameters, i.e. A* = nb* which is the characteristic nuclear 
quadrupole interaction strength of TLS's consisting of n atoms and the tunneling matrix 
element A between the left and right wells. When A > A*, the nuclear quadrupole 
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interaction can be neglected. At smaller tunneling amplitude A < A* the effective tunneling 
coupling is reduced exponentially due to the small overlap between different nuclear spin 
eigenstates in different wells. This is similar to the polaron effect and we use this similarity 
in order to suggest a solvable toy model based on the replacement of the nuclear spins with 
oscillators. 

The resonant dielectric constant is determined by the contribution of all resonant TLS's 
with a characteristic tunneling amplitude of order of energy, which varies from the thermal 
energy to some characteristic maximum value T < A < W . The deviations from the stan- 
dard tunneling model are observed at temperatures comparable with the nuclear quadrupole 
interaction, i.e., for T ~ A*. They show up as a temperature-independent plateau due to the 
breakdown of coherent tunneling in the energy range T < A < A*, where resonant TLS's 
do not exist. Therefore, TLS's with tunneling amplitudes belonging to this domain do not 
contribute to the resonant dielectric constant. Note that the hypothesis of the breakdown 
of coherent tunneling below T ~ lOmK has been proposed by Enss and Hunklinger [39] 
assuming, however, that it is due to the interaction between TLS's. Our work suggests the 
valuable realization of their hypothesys. 

The results agree qualitatively with a number of low-temperature dielectric measurements 
made by different groups. To obtain quantitative agreement for the temperature at which 
the plateau forms, we have to assume that the relevant two- level systems consist of, at least, 
four atoms (n — 4). This assumption agrees qualitatively with the TLS model based on the 
renormalization group theory [12] and is in line with molecular dynamics studies of glasses 
[22]. 

At very low temperatures T < T min (T min < O.lmK) the plateau in the resonant dielectric 
constant should go over into a logarithmic temperature dependence characterized by the 
same slope de res /d\n(T) as found at high temperatures, i.e., for T > T max ~ A*. Therefore, 
it is worth while attempting an experimental observation of this behavior. It is unclear 
whether it is possible to perform measurements at such low temperatures at present. A 
comparison of the measurements with the theory can be used to estimate the number n 
of tunneling atoms involved into single TLS because the logarithmic width of the plateau 
\ia(T max /T min ) is approximately equal n. 

The nuclear quadrupole interaction should affect the sound velocity since it enters its reso- 
nant part in the same manner as to the resonant dielectric susceptibility. Although it is clear 
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from available experiments that there is a difference in the nature of TLS's contributing to 
dielectric and acoustic properties, there are observed correlations between the two responses 
both in the hole burning and non-equilibrium dielectric measurements [3, 5, 40]. We pro- 
pose to extend acoustic measurements of glasses possessing atoms with nuclear quadrupole 
moments to lower temperatures in order to search deviations in the temperature dependence 
of a sound velocity. The most promising candidates are those materials in which the plateau 
in the dielectric constant is seen, including BK7, 10% K-Si0 2 and BaO-Al 2 3 -Si0 2 glasses 
(see Table 2). 

We have also analyzed the dependence of the resonant dielectric constant on the external 
magnetic field. We show that a high magnetic field B ~ 5 — 10T affects nuclear spins stronger 
than the nuclear quadrupole interaction and therefore will restore coherent tunneling. This 
is because the nuclear spin states in the right and left wells get aligned parallel to the 
field. We have performed the model calculations of the dielectric constant in a magnetic 
field by combining numerical simulations for the behavior of a single spin affected by both a 
quadrupole interaction and an external field. The experimental verification of our theory can 
be made measuring the dielectric constant in the strong external magnetic field, which should 
eliminate the temperature independent plateau and restore the logarithmic temperature 
dependence. 

To our knowledge, most of experimental studies of the magnetic field effect on the di- 
electric constant have been performed at relatively high temperatures, i.e., for T > lOmK 
Ref. [16]. These studies show certain similarities and also distinctions compared to our 
predictions. The major changes in the dielectric constant are observed at T ~ 5mK. The 
application of magnetic field leads to an increase in the dielectric constant in agreement with 
expectations. At a certain magnetic field the dielectric constant of the standard tunneling 
model is restored. However,the change in the dielectric constant does not show a monotonous 
dependence on the magnetic field strength. A relatively strong effect is observed at very 
small fields (for T ~ 30mK and B ~ 0.1T). So, the behavior of the system is more com- 
plicated in this temperature range. Perhaps, the relaxational contribution to the dielectric 
constant is still significant there. Then, both the phonon-stimulated relaxation [41] and the 
interaction-induced relaxation [42] of TLS's could be affected by an external magnetic field 
and our analysis would not be applicable. The experiments at lower temperatures T ~ 5mK 
can test the model proposed here and can be used to characterize the internal structure and 
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nuclear quadrupole interaction for a single TLS. 
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